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ABSTRACT 

We perform a dimensional reduction of the U(l) x SU(2) Chern-Simons bosoniza- 
tion and apply it to the t—J model, relevant for high T c superconductors. This pro- 
cedure yields a decomposition of the electron field into a product of two "semionic" 
fields, i.e. fields obeying abelian braid statistics with statistics parameter = |, 
one carrying the charge and the other the spin degrees of freedom. A mean field 
theory is then shown to reproduce correctly the large distance behaviour of the cor- 
relation functions of the ID t — J model at t >> J. This result shows that to cap- 
ture the essential physical properties of the model one needs a specific "semionic" 
form of spin-charge separation. 
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1. Introduction 



It is widely believed that the two-dimensional (2D) t — J model captures many 
essential physical properties of the Cu — O planes characterizing a large class of 
high -T c superconductors. 

The hamiltonian of the model is given by 

H = P [ E "*( C Ua + h.C.) + Jc\/-f Cl p ■ c}^C j5 ] P G , (1.1) 
<ij> 

where Q a is the annihilation operator of a spin \ fermion (in this paper called 
electron) at site i of a square lattice, corresponding in the physical system to 
a hole on the Cu site, and Pq is the Gutzwiller projection eliminating double 
occupation, modelling the strong on-site Coulomb repulsion (see e.g. [1]). 

In spite of the enormous efforts made so far, we still do not have a good mean- 
field theory for the 2D t — J model, i.e. such a successful saddle-point that we can 
apply the standard quantum many-body techniques to calculate the fluctuations 
around it for describing the essential physics. On the other hand, we have a much 
better understanding of the ID t — J model. In fact, the strictly related large U 
Hubbard model has been solved by Bethe-Ansatz [2] and for U ~ +oo the ground 
state wave function can be written [3] as a product of a Slater determinant for 
spinless fermions, describing the charge degrees of freedom, and a ground state 
wave function for a "squeezed Heisenberg chain", i.e. a chain where all empty 
sites are "squeezed out", describing the spin degrees of freedom. For a finite 
number of electrons, the pseudomomenta of the spinless fermions are related to 
the "spin rapidity" of the spin degrees freedom, but in the thermodynamic limit 
the distribution of the pseudomomenta becomes a constant [4], i.e., that of a 
free spinless fermion system. For the large U Hubbard model and for the t — J 
model at t = J, the large scale behaviour of several correlation functions have been 
computed combining Luttinger-liquid [5] and conformal field theory [6] techniques. 
Related results for t » J have also been obtained in [7], by means of a more 
standard quantum field-theory approach, using a special mean field treatment of 
the non-local, or "string" field operators. Typical contributions for an euclidean 
two-point functions are found to be of the form 

p if pn 

+ (1-2) 

(x - iv c t) a - (x + iv c t) a c (x - iv s i) a s (x + iv s t) a - 

for suitable n G Z,af,af E Q, where p is the electron density, v c and v s are 
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the charge and spin velocities, respectively. Hence the correlation functions also 
exhibit charge-spin separation. These results show that the key features of the 
ID model can be understood in terms of two low-energy excitations, the holon, 
charged but spinless and the spinon, neutral with spin |. 

Partly on the basis of analogy with the ID model, Anderson conjectured 
[8] that the physics of the 2D model can also be understood in terms of low 
energy excitations, charged but spinless (holons) and neutral with spin | (spinons). 
Depending on the statistics of the holons (sometimes called "slave-particles" ) we 
have the slave-fermion [9], slave-boson [10] and the more exotic slave-semion 
approach, advocated by Laughlin [11] (generalized to a slave-anyon approach in 
[12]). The semions are special kind of anyons (see e.g. [13]), i.e. excitations 
obeying abelian braid statistics, with statistics parameter 9 = |, i.e. the exchange 
of the field operators creating such excitations produces a phase factor ±z, instead 
of +1 or -1 characterizing bosons and fermions, respectively. 

Recently, a bosonization scheme for two-dimensional fermionic systems has 
been proposed, based on the introduction of an abelian Chern-Simons gauge field 
[14]. Such a scheme has been extended to a non-abelian version in [15] and both 
versions have been applied to the t — J model in [16] (see also [17,12]). The U(l) 
Chern-Simons bosonization has been shown to correspond essentially to the slave- 
boson and slave-fermion approaches (depending on the choice of the gauge fixing 
[16]); while the non-abelian U(l) x SU(2) Chern-Simons bosonization corresponds 
to the slave semion-approach. 

Although every bosonization scheme yields an exact identity between cor- 
relation functions of the original fermionic field and suitable bosonic correlation 
functions, the mean field approximation (MFA) gives different results in different 
bosonization schemes. It is then natural to ask which one of these schemes has a 
better chance to describe correctly in MFA the physics of the model. One expects 
that possible indications for the answer might be obtained from comparison with 
the known analytical results of the ID t — J model. 

One is then naturally led to discuss a dimensional reduction of Chern-Simons 
bosonization to ID systems. One can verify that this reduction corresponds, 
roughly speaking, to a Jordan-Wigner-like bosonization for the gauge group U(l) 
and a suitable non-abelian generalization of it for the gauge group U(l) x SU{2). 

In this paper we show that the large distance behaviour of the correlation 
functions of the ID model are indeed reproduced by a mean field theory of the 
U(l) x SU(2) bosonization; thus, interpreting spinon and holon fields as the ID 
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counterparts of semion fields, i.e. they obey abelian braid statistics with statistics 
parameter 6 = \ (in ID only the statistics of fields, but not those of excitations 
are well defined, see, e.g. ,[18].) It turns out that this dimensional reduction gives 
essentially a more systematic justification and refined structure to the approach 
followed in [7]. 

This result shows how to obtain the features of the ID model in terms of 
standard quantum field theory techniques, and encourages us to pursue the study 
of the non-abelian Chern-Simons bosonization of the 2D t — J model, suggesting 
some ideas for developing a possibly reasonable mean-field treatment. Further- 
more, it implicitly supports the interpretation of spinons and holons of the 2D 
model as semions, if they are still well defined excitations. 

The plan of the paper is the following: 

- in sect. 2 we outline the Chern-Simons bosonization scheme and apply it to the 
2D t - J model; 

- in sect. 3 we perform the dimensional reduction to the ID model for the partition 
function; 

- in sect. 4 we discuss the mean-field approximation; 

- in sect. 5 we perform the dimensional reduction of correlation functions and discuss 
their mean-field treatment. 

Some detailed computations are deferred to Appendices. 

2. The Chern-Simons bosonization 

We recall the main definitions and results of Chern-Simons bosonization 
scheme applied to spin | fermion systems in 2D [15,16]. 

Let ^q,, 1 !/* (resp. $ a ,$*),a = 1,2 be two-component Grassmann (resp. 
complex) fields describing the degrees of freedom of a spin | canonical non- 
relativistic fermion (boson) field operator if! a (resp. $ a )- 

Consider a system of spin | non-relativistic fermions interacting via an instan- 
taneous, spin independent two-body potential, and in the presence of an external, 
minimally coupled abelian gauge field A. The classical euclidean action of the 
system is denoted by S(^^*\A). Let B (resp. V) be a U(l) (resp. S*7(2)) 
gauge field and let 
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Sc.s.(B) = ^~j d 3 xe^ p B»d»BP(x) 

Sc.s.(V) = ^~j d 3 xTre^ p (v^V? + ^V^V^(x) 



(2.1) 



be the corresponding euclidean Chern-Simons actions. Then the following bosoniza- 
tion formulas can be derived: 

1) the grand-canonical partition function of the fermion system is given by: 



J VB J VQV^e-^^'** \ A + B )+ S ^s .(B)] 

/ V BVV j VQV^e-^^^* \a+b+v)+2S c . s . (b)+s c . s . (v)} 

where a) corresponds to the [/(l)-bosonization and b) corresponds to the 

£7(1) x 5'?7(2)-bosonization and gauge fixings for the respective gauge symmetries 

of the actions are understood. 

2) Let 7a;, x = (x°,x) denote a string connecting x to infinity in the ir -euclidean 
time plane; then one can prove an identity between the correlation functions of 
\I/ a , \&* in the fermionic theory and the correlation functions of the non-local fields 

a) 

$ Q ( 7x |S) = S $ a (x), K(lx\B) = e~ 1 ^ B ( x ) 

in the £/"(l)-bosonized theory, and 
b) 



QatixlB, V) = eh* B (Pe^ V ) cef3 $ p (x), 
K(l*\B, V) = ^(x)e~ l L B ( Pe - i L V ) poi 



(2.3) 



in the U(l) x S'£7(2)-bosonized theory, for spin-singlet correlation functions. P(-) 
in eq. (2.3) denotes the path-ordering, which amounts to the usual time ordering 
T(-), when "time" is used to parametrize the curve along which one integrates. 
(For a careful discussion of boundary conditions and further details, see [15,16]). 

These bosonization formulas can be derived using a Feynman-Kac represen- 
tation of the partition function and the correlation functions, expressing them in 
terms of brownian paths in R 2 . In this representation the only difference between 
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fermions and bosons are minus signs related to permutation in the order of initial 
and final points of the paths. Using the fact that the probability for two brow- 
nian paths in R 2 to intersect each other at a fixed time is zero, one can prove 
that the configurations of brownian paths appearing in the Feynman-Kac formu- 
las are braids with probability 1, and periodic b.c. in time convert these braids 
into knots. Finally, the minus signs associated with permutations and converting 
bosons into fermions are obtained from Chern-Simons expectation values of the 
exponential of the gauge fields £?, V arising from their minimal coupling to the 
matter bosonic fields, integrated over the knots formed by the brownian path con- 
figurations, following the construction of knots invariant in Chern-Simons theory 
(see e.g. [19]). 

As this sketchy explanation suggests, one can apply the same techniques to 
lattice theories, using a lattice version of Feynman-Kac formula [20,16], express- 
ing partition function and correlation functions in terms of random walks in the 
lattice Z 2 , retaining the Chern-Simons gauge fields in the continuum version, and 
provided the two-body potential contains a hard-core term. In fact, the proba- 
bility that two random walks intersect at a fixed time in Z 2 is not 0, so in order 
for the random walk configurations to be braids with probability 1, one needs a 
hard-core term forbidding intersections among random walks. 

As a result of this brief discussion, one can understand that the bosonization 
formulas can be applied to the 2D t — J model because: 

1) the Gutzwiller projection acts as a hard-core term; 

2) introducing a Hubbard-Stratonovich complex gauge field X one can rewrite the 
Heisenberg term as a standard kinetic term with minimal coupling to X (view in 
the bosonization procedure as introducing an external gauge field) plus a two-body 
spin-independent potential. 

In fact, the grand-canonical partition function of the t — J model at temper- 
ature T = hs/fl (where ks is the Boltzmann constant) and chemical potential fi 
can be rewritten [21] as: 

E t -j{p,n) = | MDI*P$D$*e- s, - j( *' r ' X ' X,) (2.4) 

with 

rP o 

S t -j(V, X, X*) = drJ2 - T X *<ij> X <i3> + [("* + 

Jo <ij> J 
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+h.C.] + + ^)**« + J2 U ^^JP*J^^ 



where the two-body potential is given by 



(2.5) 



1,3 
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+oo i = j 
otherwise. 



(2.6) 



(The euclidean-time (r) dependence of the fields here and after is not explicitly 
written, and repeated spin indices are summed over.) 

The bosonized action is obtained via substituting the time derivative by the 
covariant time derivative and the spatial lattice derivative by the covariant spatial 
lattice derivative: e.g. in the U(l) x SU(2) bosonization 



vl>* — v|> 



vl>* v[r 



d_ 



+ iB (jj)l + iV (J) 



a(3 



One can then decompose 



with the constraint 



(2.7) 



(2.8) 



(2.9) 



where E Q and E are a 2-component and a 1-component complex lattice fields, 
respectively. However, such decomposition is ambiguous, since we can still perform 
a local Z7(l)-gauge tranformation leaving $ invariant: 



E je iA i, 



(2.10) 

Sj a — > Ylj a e~ j , A, e[0,27r). 

Therefore, the theory expressed in terms of $ is equivalent to the theory in terms 
of E and E only with a gauge-fixing term of (2.10), not breaking the U(l) x SU(2) 
gauge invariance of the $-theory. As an example one can choose a Coulomb gauge 

if V 

for the U(l) x 5 , t/(2)-gauge invariant field "arg (S*Pe J <^> £_,•)" supported on 
links: 



<ij>:£e<ij> 



arg (E*Pe 



(2.11) 
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In terms of E and E the U(l) x S'L r (2)-bosonized action of the t — J model is given 
by 



^ r 2^ -jX <ij> X <i j > 

<ij> 

+ [(-t + X <ZJ> )E*e^> B E i Z* a (Pe i f<a> V ) a p^ jP + h.c] 

8 J ~ ~ ~ / 8 \ (2-12) 

j 

+ Y J Ui,jE*E i E*E j + 2S C . S .(B) + S C . 8 .(V) 

with constraint (2.9) and gauge fixings understood. 

One can easily convert E into a fermion field E, by inverting the sign of 
S C . 8 .(B); one then omits the hard-core term in Uij. (In fact to perform this 
transformation we couple E to a new U(l) gauge field B' with action S C _ S _(B'), 
changing variable B — > B" = B + B' and integrating out B the claimed result 
follows; one then rewrites B" again as B.) Integrating out X, X* and using the 
anti-commutation properties of the field E, one obtains the action 



S t -j(E,E*,X,X*,X,X*,B,V) = 

Jo 



S t _j(E, E*, E, E*, B, V) = jf drJ2 E j(^+ ^o(j) + I* + \)E 3 
+ (-tE*e^> B E i ^ ja (Pe i ^> V ) a pL ip + h.c.) 

<ij> 

+ l E *E j E*E i {m a (Pe i f«i> V ) a ^\ 2 - \) - 2S C . S .(B) + S C . S .(V). 

(2.13) 

Notice that if instead of E we use the bose field E, the J term gets an opposite 
sign due to the commutation properties of E. 

In terms of these variables the correlation functions of \1/ Q are given in the 
bosonized theory by 

d> a (7*l#, V) = eL B E x (Pe i L V ) a ^ xf3 . (2.14) 
Therefore one can view our original electron field as a product of U(l) x SU(2)- 

i f B 

gauge invariant fields E(p/ X \B) = e J ^ E x (holon) and T, a {lx\V) = 
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(Pe ) a f3^xf3 (spinon). From the coefficient of the Chern-Simons terms of B 
and V one can derive [15] that the corresponding non local holon and spinon field 
operators obey semionic statistics, i.e. their statistics parameter is 6 = \. 

3. Reduction to ID: the partition function 

The reduction of the 2D-system to ID is obtained by letting the electron field 
sit on the lattice Z instead of Z 2 . Denote by 1 the spatial dimension in R 2 along 
the ID lattice Z. 

Since in the partition function the only dependence of B2 and V2 is in the 
Chern-Simons action, one can integrate them out yielding the constraints 

i) 

e^B v = 0, 

ii) 

e liV (& i V + V' i V) = ; ^ = 0,1. (3.1) 

The Coulomb gauge-fixing condition for the U(l) gauge invariance, B\ = 0, with 
the boundary condition B^x 1 = +00) = 0, together with the constraint i) yields 

£> = o. 

The constraint ii) is solved by 

iV^x) = g\x)d li g{x), g(x) G SU{2) (3.2) 
and in terms of g and E the SU(2) gauge transformation reads as 

9j -> ''/.'// • >-., -> rjXj, S* -> E*r), (3.3) 

where rj G SU(2) and gj(r) = g(r,j). 

To gauge-fix (3.3) we use a ferromagnetic reference spin configuration, by 
setting 

=i=(J)- (3-4) 
Finally the gauge-fixing (2.11) reads, being reduced to ID: 

arg mgjgjXj) = 5. (3.5) 
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Choosing the constant 5 in (3.5) to be for later convenience, in the gauge (3.4), 
this implies that (gjgj)n is real and positive, hence gjgj is in the linear span of 
II, <r x , a y , because gjgj G SU (2). 

Let us count the degrees of freedom (d.f.): the matter fields i?,£ , *,E a ,S* 
have 2 + (4 — 1) d.f., the field g has 3 d.f. and the gauge fixings (3.5) and (3.4) 
eliminate 1+3 d.f., so we are left exactly with 4 d.f. as we had for \l/ a , ; the 
charge degrees of freedom are carried by E, E* , while the spin degrees of freedom 
are carried by g. 

The partition function of the ID t — J model can be rewritten as 

E t _j = f VgVEVE*e- s *-^ E > E *^ \[ 5(arg (gjgjU) (3.6) 

J <ij> 

with 

St- J = J P dr E* (§p + V + \ + (gfagsh^Ej 

3 (3.7) 
+ (-tE^E^g^ + h.c.) + ^EjEtE^Kgjg,)^ 2 - i}. 

<ij> 

[Weusedf slxR2 (g- 1 dg) 3 = 0]. 

We now find a configuration, g m (E : E*) of g minimizing the action for a fixed 
configuration of the holons described by E. The idea behind is that we can then 
treat spin fluctuations around g m {E, E*) in some MFA. To find g m (E, E*) we use 
the "Feynman-Kac" random-walk representation for the -E-path integration. It 
can be considered as a lattice version [16, 20] of the representation of a Feynman 
path-integration over non-relativistic fermion fields in terms of path-integrals over 
trajectories of a fixed number of Fermi particles summed over all possible number 
of particles (see e.g. [22]). 

More precisely, let A denote the ID lattice laplacian defined on a scalar lattice 
field / by 

(A/), = f i+1 + /,_! - 2/,, 
let dn(u) denote the measure on the random walks uj on the ID lattice such that 

(e~ 0A h= f (Q)=i M«>), 0>O, (3-8) 
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let Pn be the group of permutations of N elements and for it e Pn let a(ix) denote 
the order of the permutation. Then the partition function of the .E-system in a 
fixed g configuration can be written as: 

Z(g) = J VEVE*e- St - j(E ' E *' 9) = 

N=0 kEPn 3i,-,3n «j r (P)=j w{r) r=l 
N 

[II II *G&i)ii<*P { / „ drO/t^u}. 

( N J f P 1 11 

exp xzZo dT5 \^r{T)-^{T)\A\{g\^r{T))g{^3{ T ))hA 2 - d/J' ( 3 - 9 ) 

where w^ 1 (resp. uj) denotes the component of w r perpendicular (resp. parallel) 
to the time axis. We first notice that, due to the Pauli principle in ID, only the 
trivial permutation contributes to (3.9), since the random walks cannot intersect 
each other. 

Consider a fixed configuration of random walks u = {uji, ...un}', using the 
inequality |(^|^j)n| < 1 and Re(g^d T g)n = 0, one can bound from above the 
absolute value of its contribution to (3.9) in square bracket by 



N N 



n n t ex p {- J2 { 2 / dTs Mr)- Ui (r)\,i ■ i-^}- 

r=l <ij>e U ±- r,j=l J0 

This upper bound is exactly saturated by a configuration of g satisfying 

(0^0)11 = on J}, (3.10) 

(gjgj)ii = e l5 <^> , 6 <tj> = -5<^>, if < ij >e u^, (3.11) 

{g\u r {r)g{uJ 3 {T))^ = Q if \oj r (r) - ^(r)| = 1. (3.12) 
(The arbitrariness in 5<ij>, independent on time, follows from 
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by periodicity in the time direction.) 

If we further want to satisfy on uj^r the constraint 

arg (^Oii = 0, (3.13) 

we must set 5 < ij > = 0. 

To sum up, we obtain the following minimizing configuration g m = g m (uJ_): we 
choose g m constant during the period when the particle does not jump, to satisfy 
(3.10); while in a link < ij > using the representation 

g^gf = cos 91 + isen9a ■ n, (3.14) 

we see that (3.12) is satisfied by choosing 9 = ^, n z = on links not in u , whereas 
for links in ur 1 we choose 9 = to satisfy (3.11). Therefore, we can represent the 
minimizing configuration as 

gf{r\u) = e ifax ^<i^=i S " r(T) ' e . (3.15) 

(At 6 = ^, different choices of n x ,n y , but satisfying the condition + riy = 1, 
give the same element of SU(2) in (3.14); we choose to work with n x = 1, n y = 0.) 

Reexpressed in terms of the fields E, E* and denoted by g m (E, E*) the min- 
imizing configuration (3.15) becomes 

gf{r\E, E*) = e^ x e * Ei{t) . (3.16) 

(We used here the left-continuity of the paths to [20], so that at a jumping time 
uj r (r) = lim u) r (r + e), and the sign lim means e — > from the positive side.) 
It is natural to perform first the integration over g by changing variable from the 
SU ^2)-valued field g to an SI/" (2) -valued field U describing fluctuations around 
the minimizing configuration g m and defined by 

9 = Ug m . (3.17) 

For j G Z we set 

[j] = \l J ° dd ' (3-18) 

11 [2 j even, v ; 

and 
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N 

j(T\u) = J2J2 6 »AT),e- (3-19) 

£<j r=l 

To simplify the notation with the meaning being clear from the context, we set 
j( T l^) = i( r )? so that, e.g., the r.h.s. of eq.(3.15) reads e l ^ <Jx3 ^ T > . Then, in terms 
of U and the random walks u, the partition function is given by 



N=0 h,-,3N u> r (f3)=j r r=l 

J] sUgiufa)^ J! Kutu^exp f dT(Uid T U) [ar{T)]laAr) y 



<ij>E<-o r <ij>&u> r 

v 7 rP 



( . \ J f p 1 

exp ("E o dr 5 l«r(T)-w,-(T)|,i[l( c; ' t ( a; r(T))?7(a;j(r))) [5 5 r(T)][5>a . (T)] | 2 - -] J. 

(3.20) 

Analogously, defining 

j(r|£?, J E*) = X;^W, (3-21) 
we have in terms of E, E*, U: 

~ t _j = j VEVE* j VUe s t -Au,E,E*) "Q 5(arg <jj}Ui) m (3.22) 

J J <ij> 

with 

S t _j(U,E,E*) = jy T J2E*(-^ + »+^ + (U^d r U) m )E j+ 
£ [{-tE^ifJp^ + h.c] + -E^E^E^U,)^ 2 - -}. (3.23) 

<ij> 

Notice that up to now no approximations have been made and (3.20), (3.22-23) 
are exact rewritings of the partition functions of the ID t — J model. 

4. Mean field approximation 

We now wish to compare our result (3.20) with the Bethe-Ansatz ground 
state wave function of the U ~ +oo Hubbard model, essentially equivalent to the 
t — J model at J ~ +0. To make the comparison we first restrict ourselves to 
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T ~ 0, finite volume V and finite number of electrons N, with V, N being large 
and we assume p = AfyV = 1 — 5, 5 << 1 and t » J. 

Then we carry out the following mean-field treatment: 

1) We assume that the spin fluctuations can be treated in mean field in the hopping 
term of the charged particles and we denote by tn the renormalized hopping; 

2) Since the motion of the charged particles is much faster than the spin motion, 
we replace u) r (r) by its average in time, and since the paths cannot overlap 

< uj r (r) >= jz^- So we replace for the spin motion the original chain by a 
"squeezed chain" of lattice spacing (1 — and accordingly we replace 
8\u r (T)-uj(T)\,i by its mean value (1 — 5)5\ r _j\ j i in the "squeezed chain". The 
corresponding renormalized spin coupling constant is denoted by Jr. 

After making these approximations the canonical partition function decouples 
into a product of the partition function for a free charged holon system on the 
original lattice and the partition function for a spin | quantum Heisenberg chain 
on the "squeezed lattice" with lattice spacing (1 — <5) _1 : 



Z M . F .(N) = J VEVE'e-f^W^-^^**™*'^ 

5(£ E : E i - N ) ■ f VU5^g{U}U 3 ) [m ) 
i J 

where J s means that the variables integrated over belong to the "squeezed lattice" . 
(The restriction to the lattice of finite volume V is understood.) Therefore, this 
MFA correctly reproduces the features of the Bethe-Ansatz ground state wave 
function outlined in the Introduction. 

After making these MFAs one can then take again the thermodynamic limit 
of the system with a fixed density of holes. Correlation functions of fields will be 
discussed in the next section in this limit. 

To verify that the {/-system in (4.1) is actually the quantum Heisenberg chain 
one can use the CP 1 representation of U: 

where b a is a complex 2-component field constrained by 

b* a b ja = 1. (4.3) 



14 



In fact, setting 

the tZ-action becomes 

S(b, b*) = I b* lce d T b lce + 4^ ^ tfaVapkp • 6* 7 a 7<5 6 J(5 . (4.5) 

i <ij> 

One recognizes in (4.5) the action of the quantum Heisenberg chain in the Schwinger 
boson representation. The gauge-fixing condition for U becomes then 

arg b*iJ>ja = 0, (4.6) 
and it can be seen as a gauge-fixing condition for the L r (l)-gauge transformation 

b,n >b ja e«* , 6J a -6J a e-^ , <*e[0,27r). (4.7) 

(Notice that the invariance of the term with time derivative is guaranteed by the 
periodic b.c. in time and the constraint (4.6).) 

One can go a step further by using the identity a a p ■ a 7 s = — ^a^aS 

and then treating the quartic b term in the Gor'kov decoupling approximation. 
Assuming translational invariance and taking into account the 6-gauge fixing, this 
yields 

S MF (b*,~b) = ( drJ2~bj a dr~b Jct + ^ J R (b* a b ia + h.c), (4.8) 

3 <ij> 

where J R =< b*Jb ja > J R . 

In the Appendix A we show that this system is equivalent to a system of spin 
\ free fermions, described by two component Grassmann fields f a ,f*, Gutzwiller 
projected and the large-scale properties of the quantum Heisenberg chain are 
exactly reproduced by applying the standard abelian bosonization to this mean- 
field theory. The bosonization is performed in terms of a real scalar field 0_ with 
(euclidean) Thirring-Luttinger action 

S(<f>-) = ^J l(do<P-) 2 + v 2 s (d 1( P-)} (4.9) 

and with the following bosonization formulas for the left and right movers of the 
Gutzwiller projected fermions: 
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/*(*) =fi R (x) ~ (27r)-iD_(x, i) : e -**-<*> :, 

/ 2 V) =/i*(aO ~ (27r)-iD_(x,-i) : e**-<*> :, 

2 (4-10) 
=/ 2 * L (z) ~ (27r)-^_(o:, -) : e i*-(*> :, 

f 2 L (x) =n L {x) ~ (27r)-iD_(x, -i) : e-4*<*> :, 

where u s ~ Jr is the spin velocity and D_(x, ±|) is a disorder field, see [23] and 
Appendix A. 

An interesting consequence of eq.(4.10) and (A. 11) is that the spin ^ Gutzwiller 
projected fermion operator / reconstructed from the (euclidean) field / does not 
obey fermionic commutation relations but rather "semionic" with statistics param- 
eter 6 = \ . In fact, let x± = (e, x 1 ), x 1 ^ 0, let JF + (JF_) be a polynomial of expo- 
nentials of / and disorder fields with support in x° > 5 > e (x° < — 5 < — e, resp.), 
then (see Appendix A) 



]im{r-f{ t (0)f{ t (x e ±)F+) = lime-^ [arg (*±)+ ar s (-*±)l e i[«8 (-D+^g (* ± e )]. 
(r-f?(P)f?(*±)F+) = e T ^e ±|37r lnn(^_/ B (0)/ fl (xi e )J=- + > 

(4.11) 

and similar results for f L . A standard result (see e.g. [24]) of axiomatic quantum 
field theory then gives the equal-time commutation relations 

/V)/V) = e ± */(y 1 )/> 1 ), ^ * I/ 1 - (4-12) 



5. Reduction to ID: correlation functions 

In this section we discuss the dimensional reduction to ID for the correlation 
functions of the £7(1) x 577(2) bosonized t — J model. In the MFA discussed in sect. 
4 one obtains exactly the large-scale behaviour of the correlation functions of the 
ID t — J-model as derived in [5,6] by means of Luttinger liquid and conformal field- 
theory techniques. Furthermore, a simple interpretation of these results emerges 
in terms of two elementary excitations, the charged spinless holon and the spin 
| neutral spinon: the electron field operator can be decomposed into a product 
of non local holon and spinon field operator, both obeying abelian braid statistics 
with statistics parameter 9 = \. Therefore, it is natural to view (althought not 
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compulsory due to the ambiguity of excitation statistics appearing in ID) both 
ho Ions and spinons as ID analogues of semions, or, using a language more accepted 
in ID (Zamolodchikov)-parafermions of order 4 [25]. 

We start by noticing that one should choose carefully the curve 7 X , needed in 
(2.14) to define the bosonized electron field, in order to have a good dimensional 
reduction of correlation functions. In fact as mentioned in sect. 2, Chern-Simons 
bosonization is well defined only if there are no intersections in the paths on 
which the Chern-Simons gauge fields are integrated over. This was automatically 
ensured by the Gutzwiller projection for the random walks representing the virtual 
worldlines of particles, appearing in the partition function, but if we choose e.g. 
7 X as a straight line in the 1-direction 7 X may intersect many of these worldlines 
making the bosonization procedure ill-defined. This can be avoided by choosing 
7 X as the path in R 3 given by the union of the straight line going from x = 
(x°, x 1 , x 2 = 0) to x e = (x°, x 1 , x 2 = e) and the straight line joining x e to — oo in 
the 1-direction. One then takes the limit e \ 0. 

We discuss a general bosonization formula for fermion fields \l/ a , , but one 
should keep in mind that our formulas apply only if the fermion indices appearing 
in the expectation values are saturated in a spin-singlet combination, since the 
non-abelian Chern-Simons bosonization formulas have been proved only for spin- 
singlet correlations [16]. Of course, we may use the global SU(2) invariance of the 
t — J model action to reduce more general correlations to linear combinations of 
spin-singlet correlations, e.g., (** ^(x^yi^y )} = ^{^ xla (x°)^ y i a (y )). 

With a spin-singlet arrangement of indices understood and assuming x® < 
< Vs < Vs+n r ' s = 1> •••) n ~~ 1> the non-vanishing 2n-point fermion correla- 
tion functions are given in terms of E, B, E, V, by 

n n 
r=l s=l 

= (]} exp H / 5}[E x ,(^)(Pexp {-i / V})]^ (5.1) 

n 

n^»l(y;)exp{i / S}[(Pexp{z I V})X yl (y° s )}p s ) . 

S=l J ~1Vs J 1Vs 

The dependence of B2 and V2 is now both in the Chern-Simons action and in 
the strings {"f Xr ,'~fy a } attached to the fields, in the part involving the infinitesimal 
excursion in the 2-direction, as specified above. 

Let us first discuss the £? 2 -mtegration which is simpler. It yields the constraint 
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e»„d»B»(z) = Trxio.ej^C^-vi^-^-E^ 1 -^)^ -^)} (5-2) 

s r 

solved, with the Coulomb gauge fixing Bi = and b.c. B^x 1 = +00) = by 



s r 

where X[o,6] denotes the characteristic function of the interval [a, b] and © is the 
Heaviside step function. Evaluated with a regularization procedure, the term 
containing B in the action is given, in the presence of the field insertions, by 

-E'f E w(^)+E*f E ^^°) ( 5 - 4 ) 

and this is the only contribution in the correlation functions due to the field B. 
(The appearence of the factor \ w.r.t. (5.3) is due to the fact that the lattice 
points i giving non-vanishing contributions are at the boundary of the support of 
So, hence a regularization of the 5-functions involved contribute \ to the integral.) 

Remark 

Since B has its curvature concentrated at x = {x r } and y = {y s }, one can 
view the contribution (5.4) in the action as the effect of a disorder field, in the 
spirit of refs. [23,26,27]. 

Let us sketch the result of the V^-integration, and for more details see, e.g. 
[28]. To every field insertion we assign a spin | representation of 577(2) with right 
action for the creation and left action for the annihilation fields, while the Lie- 
algebra generators acting on the p-th representation are denoted by a = x,y, z. 
Integrating out V2 one obtains a matrix constraint analogous to (5.2): 



(d»V» + V»V») a {z) 



2-x[o,^ 2 ){E^ - zl ) s (y°s - ^yE*^ - - *°) (5 ' 5) 



A particular matrix valued solution of (5.5) is given by 
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(Vo) a (z) = 2n X[0 ^z 2 ){J2S(y° s - z°)Q(yl - z 1 )^ -^(x°r ~ - 



r 



{V 1 ) a {z) = 0. 

(5.6) 

The general solution of (5.5) is obtained from by an SU(2) gauge transforma- 
tion. If we perform the change of variable (3.17), it is given (at z 2 = 0) by: 



V (z) = (exp{-i^a x E$E t {z Q )}U\z))[{V ) a {z)a a + do\{u{z)- 

exp^a^i?;^ )}), 

V^z) =(exp {-i K -a x J2 E* e E e (z )}U\z))d 1 (lJ(z)e W {i\a x £ E* e E e (z )}). 

Kz 1 ' Kz 1 

(5.7) 

To understand the effect of the field (5.7), we use a random walk representa- 
tion (see [16,20] for details, with an erratum in [23]) for the correlation functions 
(5.1) at f3 ~ +oo, analogous to the one appearing in (3.20) for the partition func- 
tion: 



( ft E h(4) II E vM) U e I * B (Pe~ lf ^ V )ia r l[eLv B {Pe^ys V ) Psl ) 

r—l s=l r s 

n 



nep n k=i 4=0,1- <Mi£ (B )+<fc0)=i/i( S ) k=1 

E^i^ E Lo )= u II / vu n {^(ulu^y 

Wfa)®®} II ex p{/ dr(UHd T + V )U) MT)mT)] } ■ 

r J f p r 

exp |- Yl 2 dr5 k(T)-w'(T)|a[l(^U w ( r ))^( w '( r ))[a(T)][a'(T)]| 2 - ; 

/JO 

\ / 3 S 1 



(5.8) 
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where 

w N = {u>i, U! n , UJi,..., U>n} 

and for j G Z 

if V 

[The contribution of a path-ordered exponential Pe J ^ has been computed by 
splitting 7 2 into intervals between two consecutive crossing of (the projection in 
the — 1 plane of) 7 2 with ur. Let {zj^Zl denote the set of the spatial coordinate 
of the crossing points ordered from — oo to zi, and set zo = — oo, z p = z. Then in 
(5.8) we find: 



j=0 j=0 3 

p-1 

U (U z] (z°)e^ Ul (z )) U z i (z°)e^ s1 ^ . ] 

Apart from the presence of V and the exponentials due to the 7-strings, the key 
difference in (5.8) w.r.t. (3.20) is the appearence of n new random walks Uk, 
starting at times x® at the points x\ and ending in the points at time y®( a y, 
after wrapping = 0, 1... times around the circle of lenght (3 in the time direction; 
these paths describe the virtual worldlines of the charged particles created and 
annihilated by the insertions of the E and E* fields. Using the techniques of 
[19,28], one can show that in the presence of V every crossing of ujjj^ with 7x r (7y s ) 
contributes to (5.8) a factor ^i. 

Collecting all the factors due to the field V (5.7), we derive that 

+ iUj(r)a x U^T) (5.9) 

is the contribution of V due to the intersection at site j and time r of the curves 
lx r (ly a ) with the virtual worldlines uj_ n of the charged particles described by the 
field E. 

Combining together (5.4), (5.9) and (3.16) we find that, for the spin singlet 
correlations, the fermion fields can be exactly represented in terms of U, E, E*, as 
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(5.10) 

To obtain a more tractable expression to compare with the results obtained by 
Bethe-Ansatz and Luttinger liquid techniques, one applies the MFA discussed in 
section 4. 

An important result of the MFA is that the spin degrees of freedom appear 
in the "squeezed Heisenberg chain", where the Gutzwiller projection can be im- 
plemented exactly as a single occupancy constraint. The spin fluctuations b a in 
the squeezed chain can be then converted into fermion fields f a by 

b ja (r) = e i ^^ f ^ {T) f ja (r). (5.11) 
From the Gutzwiller projection we also derive, see (4.14), (4.25): 

/;, <' ;:, /,2. /* 2 = e^7,i. (5.12) 
Denoting quantities in the squeezed chain by [ • ] s , we obtain in MFA 



Ply 1 ] MFA 



( e ^E e<v i Ut* x U}(y°)E* e E e (y°) \ 



(5.13) 



To derive the large-scale properties we apply to (5.13) the abelian bosoniza- 
tion, obtaining at long distance, apart from an overall ultraviolet renormalization 
(see Appendix A): 

V y Vy J )l{2}yi MFA J1 V ' 

(In (5.14) for brevity we introduced the following notation: if strings of the form 
a{b} with mathematical symbols a, b appearing on both sides of an equation, the 
meaning is that the equation is valid either if we use everywhere the symbols before 
{•}, or if we use everywhere the symbols inside {■}.) 
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To extract the large scale properties of the charge degrees of freedom we apply 
the abelian bosonization also to the field E: first we rewrite it (and its conjugate) 
in terms of left and right movers: 



Ej = e ™0-&)3Ef + e'^-^Ef, (5.15) 

then to the corresponding continuum fields E L (x), E R (x) we apply the bosoniza- 
tion scheme rewriting their correlation functions in terms of a real scalar field, <p c , 
with (euclidean) Thirring-Luttinger action 



J d 2 x[{d^ 



Wc) = ^ / d 2 x[(d Q( f> c y + vi(d 1( f> c y] (5.16) 



and with the identifications 

E R = D c (x,l) :ei^ x) : etc., 

where v c ~ tn denotes the charge velocity. 

Evaluation of the E— string in the scaling limit gives 

e ±±<t> c (x)±i% (l-S)x 1 ^ (5.17) 



By the same arguments as those used in (4.11-12), one can show that the non- 
local field operator reconstructed from E(x)e ±l ^' -°° E E +E E '^ dy obeys an 
abelian braid statistics with statistics parameter 9 = \. 

Combining together (5.12-17) we have shown that in MFA the large scale 
behaviour of the spin-singlet correlation functions can be derived using the iden- 
tities: 

*; il{2} (x°) M ^( e " 1{ " 3}lf {1 - S)xl E* L (x) + e 3 ^^ xl E* R (x)y 

j D c ( X5 _l)( e l{3}if(l-<5)x 1 . e 3{l}|<Mx) • +e 3{l}*f (l-5)x 1 . e _ 3 {-l}^ c (x) .y 

D-(x,-{+}±):ei+-M: 
VH2}(V°) ~ fe 1 < 3 > 4 *< 1 -*)» 1 S L (y) + e- 3 <- 1 > < *^»' 1 £; a (y)). 

L ' MFA \ J 
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"{ + }*§ :E* R E R +E* L E L :(z)dz .R{*L\ , N 

£> c (y,+l)(e 1{3} ^ (1 ~ <5)yl : e ~ 1{_3} *^(w) : +e -3{-i}*f (i--*)?/ 1 . e 3{i}^ c (v) .). 

2?_( y , + {-}I):e-^-(«): (5.18) 

According to the general results of £7(1) x SU(2) Chern-Simons bosonization, (see 
(3.13)), the original fermion field can be decomposed into two non-local semion 
fields, £7(1) x 5£7(2)-gauge invariant: E{^ X \B), E a (7 x |y) which one may call holon 
and spinon field, respectively. Equation (5.18) proves that for large scales, in the 

% [ Xl :E* H E R +E* L E L : 

MFA one can identify E(^ X \B) as the lD-semion field E(x)e 

and ^\{2}{lx\V) as the ID chiral semion field fi* , a Gutzwiller projected chiral 

fermion field of the squeezed Heisenberg chain. 

In Appendix B we show that the formulas (5.18) when applied to the corre- 
lation functions of the ID t — J model in the regime t » T, indeed reproduce 
correctly large scale behaviours identical to those obtained with Bethe-Ansatz and 
Luttinger liquid techniques extrapolated from the large £7 Hubbard model [5] and 
the t = J, t — J model [6]. 

To conclude, we have shown that one can obtain the correct large scale be- 
havior of the ID t — J model by simply using a mean field theory treatment of the 
dimensional reduction of the U(l) x 577(2) Chern-Simons bosonization. Moreover, 
we have shown that the U(l) x 577(2) Chern-Simons bosonization is the most nat- 
ural mathematical framework for describing the spin-charge decomposition of the 
electron field in terms of semionic fields. This shows that the key physical proper- 
ties of the ID t — J model are captured not by an arbitrary spin-charge separation 
scheme (in fact slave fermion and slave boson approaches failed to reproduce the 
correct large scale behaviour of the correlation functions) , but rather by a specific 
semionic form of the spin-charge separation. This gives rise to our hope that the 
power of this formalism and the underlying physical intuition will survive in 2D. 

Appendix A 

We first fermionize the system (4.8). Since the spin components of b a are 
coupled via the Gutzwiller constraint, one cannot apply to them independent 
Jordan-Wigner transformations. To derive the correct transformation we use once 
more the reduction from a 2D system. Coupling the b field to a £7(1) gauge field B 
with action S C _ S .(B) in 2D is known to convert b to a fermion field /. Integrating 
out B we get the constraint 
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e ltv & t B v (x) = 2nY J Khc t {x°)5(x 1 - j), = 1,2. 

j 

In the Landau gauge, d^B^ = 0, it can be solved by 

B A X ) = ^^*a^a<> )^arg (x 1 - j). 

j 

Reducing this formula to ID we find 

B (x) = 0, B 1 (x) =7Tj2~b* j aha(x°)d i e(x 1 -j). 

j 

Applying the Gutzwiller constraint we obtain 

B 1 (x)=nj2^ 1 -J), 

j 

if B 

so that e J < i i> = e l7r and the action for fermions reads as 

s(f, n = [drj2 f* d rfi - E Mfifi + h - c -) 

i <ij> 

with the constraint 

flafja = 1, (A.2) 

i.e. it describes exactly a system of spin \ non-relativistic free fermions, Gutzwiller 
projected. 

It is well known [29] that the large-scale properties of the spin \ quantum 
Heisenberg chain can be described in terms of a real scalar field ip with Luttinger- 
Thirring action 

^M = £ J[(d 0V ) 2 + v 2 s (d lV n (A3) 

where £ = 2 and v s denotes the spin velocity. (Here we use the convention in 
which a mass term would be described by : cos <p :.) 

Let us now show that we indeed recover this result starting from the mean 
field model (A.l). To analyze the large scale properties, following [30], we first 
introduce the decomposition of / in right and left movers in a lattice labelled by 
sites in (2Z + \): 
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f2na -I n /( 2n+ l) Q + ( l) n /(2n+i)a' 

f _,-2n+l f L , / -\2n+l f i? 

J2n+la — I J(2n+±)a^~y L ) J {2n+\)a' 

The Gutzwiller constraints are given by: 

/*L rL I ("*_R eR -i 

j-*_R j-L , r*L eR — ' 

J2n+±aJ2n+±a •> 2n+\aJ 2n+±a ~ U- 

The fields / L , f R are then assumed to have a good continuum limit with linearized 
dispersion relations, resulting in a large-scale continuum action given, before the 
constraint is implemented, by 

s MF (f,n = J d 2 x[f: R (do+iv s d 1 )f R + f: L (d -iv s d 1 )f^. (ag) 

Introducing two real scalar field 4> a ,cx = 1,2 and applying the standard 
abelian ID bosonization (proved in [23,31] to be a special version of the dual- 
ity transformation), we obtain a Luttinger-Thirring action 

2 

s MF ((p a ) = J2 s ^> ( A7 ) 

a=l 

and the bosonization formulas for fields [23,27]: 



(A8) 



f*{x)~{2it)--*D a {x,l):e-i^ :, 
f* R (x) ~(27r)-^D a (x,-l) : ei*^ :, 

/ Q L (x)~(27r)-^ Q (x,l):e^«^ :, 

r a L {x)~{2n)-^D a {x,-l):e-^^ :, 

where D a (x, ±1) is a disorder field and we adopted the notations of reference [23], 
to which we refer for more details. [In (A.8) : :,(eR denotes the normal 

ordered exponential, defined as follows: let 5% be a regularization of the Dirac 
5-function S x , with u.v. regulator k, i.e. k^oo^x] then 



with 



: e «*<*> := lim e ^(^) (27r) C 2 e -^(^,A-^) {A g) 



A 1 (x, y) = — In yj{x° - y ) 2 + v^x 1 - y 1 ) 2 . 
2tv 



Formally, we rewrite (A. 9) as 
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. e <j>{x) := e^(-)(27r)^e-^ A " 1(a; ' a;) . ] (A10) 

In the model with action (A. 3) the expectation values of products of disorder fields 
and exponentials are given by 



< 



l[D(x J Xj)ef' l,ix)fiv)dv > 



0, if Ej0^0or/(0)^0 for /real, 

e IE i </i0 ln |a ; i-^l e / d 2 xd 2 yf(x)ln\x-y\f(y)^J2 j J d 2 xf (x)a,vg [(x° -x°)+iv s (x 1 -x))] 

(All) ! 



where \z — w\ = a/ (z° — w ) 2 + v 2 (z x — w 1 ) 2 . 

One can implement the Gutzwiller constraint exactly, defining 

01 ± 02 

Then, with the notations of (A. 3 ), we have 

Smf(<P+, <P-) = S 2 (<p+) + S 2 (<f>-) (A12) 

and, e.g., 

f?{x) = D + (x, !)£>_(*, 1) : e-i^+W : e"^"^ : 

etc, so that (taking care of the normal ordering, see [27]) the constraints (A. 5) 
become 

^- ■■di<P+ := 0, 

271 n (A.13) 

: cos (cp + H ) :: cos 0_ := 0, 

2 

solved by 

<f) + = const, : e i<t>+ := 1. (A14) 

As a conseguence Smf{4>) = ^ (</>-), i-e. we recover eq. (4.9) and we obtain the 
bosonization formulas (4.10). 

From equations (4.10) and (A. 11) we derive, e.g., 
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(/f(x)/ T (y)> = e^ 1 -y 1 Hf* R (x)f 1 R (y)) + e-^ 1 -y 1 Hf^(x)f 1 L (y)) = 

(x 1 - y 1 ) - iv s (x° - y°) ^{x l - y 1 ) + iv s (x° - y°) 

Therefore the mean field treatment of Gutzwiller projected free electrons repro- 
duces exactly the large scale properties of the Heisenberg chain. 
Remark : Suppose we keep the hard-core constraint for the individual fields b a , 
but perform a mean field treatment of the remaining Gutzwiller constraint, then, 
since the two components of b a are not any more coupled we can fermionize them 
by separate Jordan-Wigner transformations. To the corresponding fermionic field 
f a one can apply the same treatment as before, but since the Gutzwiller constraint 
disappeared, our bosonization formulas are simply eqs.(4.18). This procedure re- 
produces the approximate results of [7] . Therefore, the treatment of the constraint 
should be exact, and only afterwards one can use MFA. 

Finally, let us apply equations (4.10) and (A. 9) to prove (5.14): we rewrite 
(5.13) as 

e -(-)"i4«(i-*) e 4* + (») £>+(y , |)- 

(V*§ (l-^ 1 . e -i<t>+(y) .. e (-fi<P-(y) . +e +i§(l-S)y 1 e ^^ + (y) .. e -(-)"i*-(tf) 

(A15) 

with the Gutzwiller projection implemented imposing equation (A. 14). Formally, 
see (A.9-A.10), we have 

e i<t>+(y) :e -i<t>+(y) >e 4 A "'(w) ) 

e i<t>+(y) :e i4>+(y) gfA- 1 ^^) . e i4>+(y) e |A- 1 ( y , y )_ 

Therefore, the ratio of the coefficients of the second to the first term in (A. 15) 
tends to as the U.V. cutoff (assumed to be the inverse of the scale parameter) is 
removed. Recalling eq.(4.10), at large scale we recover (5.14). 

Appendix B 

Using (5.18) and Appendix A we compute the large scale behaviour of 
1) density-density correlation function: 

(: vJr Q : (x) : : (y)> ~ <: E* E : (x) : E*E : (y)> 
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1 1 

~ : di<t>c ■■ (x)— : d l( j) c : (y)) 



An 2 



+ 



I 1 

;! - v 1 ) + iv c (x° - y )) 2 J 



(( a .l_ y l)_ it , c ( a .0_ y 0))2 (( a -l_ z/ l) + ^ c ( a -0_ z/ 0)) 
27T 2 (X 1 -^1)2+^2(2.0 _ y 0)2' 



2) spin-spin correlation function 



<*^M*)*;^r**(y)> = 



MFA 



l(E*E(x)f^f 1 R (x)E*E(y)nf 1 R (y))+ 1 -(E*E^^^ 
\ E *E{x)e^-~- E * RER ^ 



2 

+ 



4 2 

(ri R (x)fi R (y))(f*i L (x)fi(y)) + {e -™fl!-- E * RER + E * LEL -^ d *). 

e- i ^ 1 - S){yl - xl) (f R (x)r R (y))(f L (x)r L (y)) 

= <j2 cos [(njl-SW-x 1 ))] 

({x 1 -y 1 ) 2 + v 2 (x° -y^dx 1 -y 1 ) 2 +v 2 (x° -y ) 2 )^ 

where in the third equality we use E*E(x) =: E*E : (x) + 8 and we neglect the 
first term as subleading; 

3) electron-electron correlation function 

(#*(:r)#a(y)) ~ (( J E* L (x)£; L (y)e-^ (1 - ,5)(xl - yl) + J E;* i? (x) J E; i? (y)e 3 ^ (1 - ,5)(xl - yl) )- 



e 



-ii/; i :^^ +B -^: Wd , /rfl(a . )/iB(y)) 

+ (( J E* L (x)£; L (y)e-^^ (1 - 5 ^ xl - yl) + ^(z^^y)^!- 5 )^ 1 -^)). 
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3 if ti-.E-«E« + E--Ei*z)dz f L ix)f *L( y y = - 1)Dc fa 1)D _( X , -\)D_{y, \) 



e -if (l-5)(x 1 -y 1 ) . e |<^ c (x) .. e - j0 c ( V ) . +e i3f (l-5)(x 1 -y 1 ) . e ~i\4> c {x) .. e i|0 c (y) . 

1 1 



)+/l.C. 



[(a; 1 — y 1 ) + zf s (x° — y )] 2 [(x 1 — y 1 ) 2 + v 2 (x° — y ) 2 ] is 

+ /i.e. (S.3) 



- + 



We see that the large scale behaviour of the correlation functions is of the form 
(1.2) and a comparison of the values of n, af, af with those found in [5] and [6] 
shows that they agree exactly with the result obtained for the large U Hubbard 
model and the t — J model at t = J, extrapolated to the region t » J. 

Acknowledgments. Useful discussions with J. Frohlich are gratefully acknowl- 
edged. 
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